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Abstract. An "exact discretization" of the Schrodinger operator is considered 
and its direct and inverse scattering problems are solved. It is shown that a 
differential-difference nonlinear evolution equation depending on two arbitrary 
constants can be solved by using this spectral transform and that for a special 
choice of the constants it can be considered an integrable discretization of the 
KdV equation at large times. An integrable difference-difference equation is also 
obtained. 
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1. Introduction 

We are interested in discrete dynamical systems that in the continuous limit give the 
KdV equation and that are integrable by means of the inverse scattering technique. 

It is long time that the inverse scattering method was extended to solve 
differential-difference and difference-difference nonlinear equations. See in particular 
the pioneer works of Flaschka [Q , Manakov Q , Ablowitz and Ladik ||, Q , Kac and van 
Moerbeke j^]. In these papers and in following papers (among other discrete equations) 
different versions of the discrete KdV equation were introduced and studied. See ||- 
H , the review paper jio) and references quoted therein. 

The various already considered discrete KdV equations were related to a reduced 
form of the Ablowitz-Ladik spectral problem or to different discretized versions of 
the Schrodinger spectral problem. See, in addition, for the spectral theory of these 
discrete Schrodinger operators JllJ]-[jl3| and ju], [l5) for the case of cellular automata. 

Here we use the "exact discretization" of the Schrodinger equation considered 
recently by Shabat 

= g n 1pn+l + Xlpn (1-1) 

and obtained iterating the Darboux transformations Jl6| . 

We completely solve the direct and inverse scattering problems in the case of a 
potential g n satisfying 

+00 

(l + |n|)|s„_i -2| < 00. (1.2) 

n— — 00 

Then we associate to the discrete spectral problem ([O]) the differential-difference 
equation (p and q are arbitrary constant) 

(g n 9n-i)t = p ( — — ) + q (gn - g n -i) g n g n -i (1.3) 

\9n+l 9n-2/ 

and we solve its initial value Cauchy problem in the standard way by using the inverse 
scattering method. 

Different possible continuous limits remain to be explored. In the special case 

p = 2q (1.4) 

g n = 2 + h 2 u{nh) (1.5) 

x = nh, T = h 3 -t (1.6) 

8 

in the limit fe-iOwe get the KdV equation 

u T = -u X xx + 6u x u. (1.7) 

Therefore the equation (|l.3]), in the case p — 2q, can be considered an integrable 
discretization at large times of the KdV equation. 

Finally, we consider also a double discrete equation, i.e. discrete both in space 
and in time, associated to the double discrete Schrodinger equation 

1pn+2,m = gn,m"4>n+l,m + A?/Vi,m, (1-8) 



if we let 



that in the limit of continuous time reduces to (1.2). 
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2. Differential-difference equation 



In the following it is sometimes convenient to use, for any discrete function /(n), the 
notation 

fi = f(n + i). (2.1) 

when n is considered generic but fixed. 

Then the discrete Schrodinger equation (1.1) reads 

ip2 = gaipi + A^o- (2.2) 

We consider the auxiliary spectral problem 

i> xt = (A + XB ) + (Ci + ADi)Vi, (2.3) 

and we impose compatibility in order to determine the coefficients in the r.h.s. and 
the nonlinear differential-difference evolution equation which can be solved by using 
the spectral transform related to (Im- 
precisely, let us denote by E the shift operator, Ev n — v„+i- Then, the 



compatibility condition between fl2.2| ) and (2.3), that is Eipu — dt ip2, gives 

(ill +\B 1 )i> 1 + \ (C 2 + \D 2 ) + (C 2 + \D 2 ) g A 

= (A-j. + AB_i) Vi - (i4_i + AB_i) 5 _i^ + A(C + A£» )^o 

+ g ot ^i + (A + XB a )g ^ a + (d + XD^g^. 

Assuming ip\ and tpo to be linearly independent and equating the coefficients of the 
different powers in A we get the following relations 

D 2 = Do (2.4) 

C 2 -C =g B -g- l B_ l (2.5) 

A Q g = (2.6) 

B 1 - £_! = g {D x - D 2 ) (2.7) 

ff ot =A 1 -A- X +g {C 2 -C 1 ). (2.8) 



From (2.4), 2.6) and ( |2.7| ) we immediately have 
D(rc) = rf 
A(n)g(n) = p 
B(n) = q 

where p, q and d are constants, and then we are left with 
C 2 - Cq = q (g Q - g_i) 

and 

1 1 

9-1, 



got =p 



+ go(C 2 -d). 



We write ( 2.13 ) for the shifted indices 



g-u = v ( — - — ) + 9-i{Ci - Co) 

,5o 5-2/ 



(2.9) 
(2.10) 
(2.11) 

(2.12) 
(2.13) 

(2.14) 



and sum it to ( p. 13 ). Using (2.12) we finally get the searched differential-difference 
equation 



(gog-i)t = p ( — - — )+ q (go - 5-1) .905-1- 

51 5-2/ 



(2.15) 
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of the ai 
A(n) = 



The coefficients of the auxiliary spectral problem are given by 

P 



g( n ) 

Bin) = q 



(2.16) 
(2.17) 
(2.18) 



C{n) =qJ2(9(n - 2j) - 9(n - 2j - 1)) 

D(n)=Q, (2.19) 

where for convenience we chose equal to zero the integration constant of C and D. In 
the following we need the limit of these coefficients as n — ► ±oo. For n — > — oo we 
have C{n) — > 0, but for n — » +oo it is necessary to consider separately the cases odd 
n and even n. We have 

m— 1 

C(2m) = g ^ (ff(2r) - fl (2r - 1)) 



C(2m+l) = g ^ (ff(2r + l)- fl (2r)). 



The limits 



lim C(2to) 

m — > + oo 

lim C(2m + 1), 

m — ^+oo 



thanks to (1.2) exist and satisfy 
c e + c d = 0. 

3. Spectral transform 

3.1. Direct problem 



(2.20) 
(2.21) 

(2.22) 
(2.23) 

(2.24) 



(3.1) 
(3.2) 

(3-3) 



The spectral problem (1.1) by using notations 

A = -1 - k 2 

g(n) = 2 + u(n) 

ijj(n;k) = (1 + ik) n x(n;k) 
can be more conveniently rewritten as 

(1 + ik) X (n + 2; k) - 2 X (n + 1; k) + (1 - ik) X (n; k) = u(n) X (n + 1; k). (3.4) 
The Jost solutions can be defined via the following discrete integral equations 

I - >k\ j ~ n 

- " 5 «C?-i)/i + O';*0 (3-5) 



/i (n; fc) = 1 + 



2ik 

1 

2ik 



j=n+l 



E 



i 



i - 



1 -ik 



1 + ik 
1-ik 



j-n 



u(j~l)n~(j;k). (3.6) 



] = -co 

We show in the next section that for a potential u(n) decaying sufficiently rapidly to 
zero at large n the Jost solutions a + and are, respectively, analytic functions of 
k in the upper and in the lower half-plane and have a continuous limit on the real k 
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axis. Then, it is natural to introduce the following spectral data (defined on the real 
axis fccj = 0) 



-. +00 



j = -oo 
+00 



6± (fc ) = ± J_ y I 

w 2ik VI 



1 + ik\ j 

— ik 



«Ci- 1)^(7; *) 



P ±(k) = 



b ± (k) 

a±(fc)' 



(3.7) 

(3.8) 
(3.9) 



Due to the above mentioned properties of analyticity of p^ , it is clear that a can be 
analytically extended, respectively, to the upper half-plane and to the lower half-plane 
of the complex spectral parameter. The function 1/a plays the role of transmission 
coefficient and p is the reflection coefficient. 

By using the definitions of /i ± in ([T^) and ( |3.6| ) one can check directly that on 
the real fc-axis, where both Jost solutions are defined, 



V + {n;k) 
a+(k) 



- p (n; k) 
1 ™ 

— y 

2ik ^ 



1 - 



1 + ik 
1 - ik, 

l + ik^ n -' 



= P + {k) 

u(j - 1) 



1 + ifc V \ p + (j;k) 
1-ik) [ a+{k) 



v 1 — ik J \\—%Kj crnft) 

j = -oo L 'J 

so that, assuming that the "integral equation" (|3.6|) is uniquely solvable, we get 



p + (n;k) 
a+(k) 

Analogously we obtain 

p~(n; k) 



p (n; k) + 



p + (n; k) + 



1-ik 
1 + ik 

1 - ik 
1 + ik 



a-{k) 

Taking into account the integral equations for p^, we see that 

lim p + (n; k) = 1, > 0, 

n — >+oo 

lim p~(n;k) = 1, /cq < 

n — > — 00 

while for fcc$ = we have 

1-ife 



p + (k)p~{n;-k), fc = % (3.10) 



/? (fc)/i + (n; — fc), fc = fcsR. (3.11) 



(3.12) 
(3.13) 



/^(n; fc) ~ a ± (fc) + 



1 + iA; 



6 ± (fc), for n -> too. (3.14) 
= 1 for k — /csr. Then, 



l-ifc 



1+ifc 



Note that ( ) gives rise to an oscillating term, since 

from the asymptotics at large n of (B.1C) and ( |3.1l| ) we get the unitarity relations 
(k = kn) 

a + {k)a-{k) + b + (k)b-{-k) = 1 (3.15) 
b+(k)b-(-k) = b+(-k)b-(k) (3.16) 
a+(k)b-(k) + b+(k)a-(-k) = 0. (3.17) 
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v + (n; k) 



v (n; k) 



1-ik 
1 + ik 
1 -ik 



fi (n;-k) 



fi + (n; -k) 



(3.18) 



(3.19) 



l + ik / 

which are solutions of the spectral problem fl3.4|) and, up to the kinematical factor 
jt^ J , are, respectively, analytic in the upper and lower half-plane of the spectral 
parameter. 

If one defines the discrete Wronskian as 

'1 + ifc N " 



1-ik 



[x(n + 1; k)(p(n; k) - x(n; k)ip{n + 1; k)} 



(3.20) 



it is easy to verify that for two solutions \ an d <p of (f^J) the Wronskian is independent 
of n. Then, we can evaluate W(fj, + using the asymptotic limits of /i + (n;/c) and 
H~(n; —k) for n — » — oo in ( |3.12j ), ( 3.13| ) and ( 3.14 ) getting 

W(n+,v+)(k) = j^a+(k). (3.21) 

Of special interest are the zeros of a + (k), say at k = fc+ (m = 1, 2, . . . , N + ), which 
are located in the upper half-plane and that we expect to be related to the soliton 



solutions of the evolution equations related to the spectral problem ( |3.4[ ). From (3.21) 
we deduce that /^ + (n;£;+) and ^(n; fc^) are proportional, i.e. 



(3.22) 



The constant £>+ can be deduced by rewriting the integral equation (3J)) for fi + at 
k = ki as follows 



l + ikj 
1 %kf] 



/i+(n;fc+) = 



T E 



1 



2ik 



1 - 



+ 00 

E 



1 — ihm, 



j = -oo 

1 ik^L \ 



1 -t~ tkm 



u(j - i)v + {j;k+) 

l + ik+Y 



u(j - 1) 



Comparing with the integral equation (3.6) defining /i we get 



D+ = 



1 



2z 



E 



l + iki 



u(j - l)v + (j;k+) 



(3.23) 



v J il'n 

Notice that the convergence of the series is guaranteed by the proportionality between 
/Lt + (ro;fc+) and v + {n]k^ n ) that have, respectively, good behaviour for n — > +oo and 
n — > —oo. 



3. 2. Neumann series for /i* 
We consider the series 



-(n;k) = ^2^ l \n;k) 



(3.24) 



z=o 
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whose coefficients are denned through the iterative algorithm 

+00 

^ +1 V;fc)= E M(n,j;k)u(j-l)^ l \j;k) 

j=n+l 

starting at fi^(n; k) — 1, where 

1 • iky '■' 



M(nJ;k) = -±- 



1 - 



1-ik 



(3.25) 



(3.26) 



If the series is uniformly convergent it solves the discrete integral equation (3.5) 
defining the Jost solution fi + . This study is most conveniently done following the 
procedure outlined in the continuous case in ]l7| ]. 

We have for j > n + 1 and fcs > 0, at the same time, 



\M{n,j;k)\ < 
\M{n,j;k)\ < 



1 



2\k\ 
1 



1 



1 + ik 



1 - ik 



j—n—l 

E 

r=0 



1 - ik 
1 + ik 



< 



1-ik 



\k\ 

< U ~ n). 



Now setting 

I(x) = 1 + \x\6(-x) 
and noticing that for (j — n) > 

U - n) < I(n)I(-j) 



we can write 



where 



\M(n,j;k)u(j-l)\<a(n)p(j) 



a(n) = ( J(n) 
0(j) = \u(j-l)\ 



1 



(3.27) 
(3.28) 
(3.29) 

(3.30) 
(3.31) 



Notice that 



a(j){3(j) = \u(j-l)\ 



(3.32) 



\k\~i 

(i + lil) ■ 

Using these bounds in ( 3.25| ) we obtain as majorant jEt® of \fjfl^\ the corresponding 
term of a series expansion of the solution of 

+00 

£(n;fc) = l+ E a(n)l3{j)Kr,k). (3.33) 

j=n+l 

One can prove by induction that for I > 

9'(„ + i,j-i) 



where 



j=n+i+l 



Q(n,j) = E a ( s )^ s )' 



(3.34) 



(3.35) 
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+ OC 



V-l 



(j + l,r-l) 



(i-l)l 
Q'-^' + l.r-l) 



(1-1)1 



j— n+1 r— j+i 

Changing the order of the sums we obtain 

-\-oc r — l 

n il+1) (n)< E a(n)p{r) £ a(j)/?(j) 

r— n+/+l j— n+1 

and then, remarking that, since r — 1) > Q(j + 1, r — 1), 
Q'0',r - 1) > Q'Cj.r- 1) - Q'(j + l,r- 1) 

- (Q(j, r - 1) - Q(j + l,r - 1)) £ Q'-^^j, r - l)Q s (j + 1, r - 1) > 

s=0 

>la(j)P(J)Q l - 1 U + hr-1) 



we derive fl3.34|) . 



Then, we change the summation limit in (3.34) as follows 



# ,+1) (n)<5>(»)/9(5)- 



/! 



and sum up over I getting 



/*(») < 1 + E "("^(j) ex P + !- i - !) 



and finally 



+ 00 

£(n)<i + c£>Cj-i)|{ 

j=n 



I(-j)I(n) ' 



where 



c = cx P £ (i + lil)l«0'-i) 

Thus for fccj > and for a potential u satisfying 



J2 (1 + \j\)\uU - l)\ <oo 



j = -oo 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



(3.40) 



the series in the r.h.s. of (3.24) is uniformly bounded and each iterated term is 
absolutely bounded by the Z-th term of a convergent series. Since the fj,"'(n; k)'s are 
analytic functions of k for fccj > 0, we conclude that fi + (n; k) is analytic for fcg > 
and has a continuous limit as kc$ — > 0. A quite similar analysis shows that n~{n; k) is 
analytic for fcs < and has a continuous limit as fccj — > 0. 
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3.3. Inverse Problem 

Let us suppose that a + has simple zeros at kj and let us denote by the residuum 
of l/a + , that is 

(k - kf) 



lim 



a+(k) 



3 = 1, 



Then we can write (3.10) in the form (k = fcjt) 



N+ 



p + (n;k) ^ + p+(n;k+) 
a+(k) 



" 3 k- k+ 

3 = 1 1 



JV+ 



M {n;-k)p + (k) 



/i (n; fc) — c 
1 - ifc 



./i+(n;fct) 



(3.41) 



. 3 k-k+ 
3=1 -J 



1 + ik, 

The terms in the left-hand side are the limiting values at the real axis of a sectionally 
holomorphic function; then, taking into account the asymptotic behavior of fi at 
large k and recalling (3.22), we can use the Cauchy-Green formula and get for fcs > 

= i + £c +M ~ (n; ~ fcj+) ( x ~ ik t\ n 

3=1 



H + (n;k) 
a+{k) 



k-k+ 



l + ikl 



1 

2i7T 



(n; — s)p + (s) / 1 — is\ 



fc-iO 



1 



(3.42) 



and for < 



JV+ 



M -(n;fc)-l + ^C+ 

3=1 



(n; — fet) / 1 - ik - 1 



1 + 



1 

2i7T 



/i (n\—s)p + (s) f 1 — is\ 



s-k + iO 



1 + is 



ds 



where 



4 D t 



(3.43) 



(3.44) 



Equations ( 3.42 ) and ( p. 43] ) allow to reconstruct the Jost solutions, given the 
spectral data. In order to complete the formulation of the inverse problem, we have 
to reconstruct the potential u. By inserting the asymptotic behaviour of p~ at large 
k 



(n;fc) = l+ ( , +o(k- 1 ) 



into (3.4) we get 



u(n) 



M(_i)(« + 2)-M(_i)(«) 



(3.45) 
(3.46) 



Therefore from (3.43) we have 

N+ 



u(n) 

3=1 



, l-ikj 



i 

2^ 



-foo 



/U (n + 2;— s) 



1 — is 
1 + is 



H (n;-kf) 



fj, (n;-s) 



1 - ifel 



1 — is 



C 



p+(s)ds. (3.47) 
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Now we want to find the time evolution of the spectral data. First, by using (3.1) and 
( |3.3| ) we rewrite the auxiliary spectral problem in terms of x as follows 

(1 + ik) xu = (A -B Q - k 2 B ) xo + (1 + ik) (Ci -D t - k 2 D,) Xl . (3.48) 

Since the /r^'s have at large n functionally different asymptotic behaviours there exists 
a x{t, n; k) satisfying ( |3.4[ ) and (|3.48| ), and an Q(t; k) such that 

X (t,n;k) = n(t;k)n + {t,n;k), (3.49) 

that is satisfying 

(1 + ik) (ntfl- 1 ^ + /i+ ) = (A -B - k 2 B Q ) (i+ 

+ {l + ik){Ci-D l -k 2 D 1 )ii+. (3.50) 

We evaluate (3.50) separately for n = 2m — > +oo and for n = 2m + 1 — > +oo using 
( |3.12[ ) and ( |2.16| )-( ^19| ) and re callin g that g(n) -> 2 as n — > oo. Then, we sum up the 
two limits taking into account (2.24) getting 

q(l - ik) 



r^rr 1 = 

2(1 + ik) 

and, consequently, the following time-dependence for f2 

P 



Q(t; k) — fio(fc) exp 



2(l + ifc) 



g(l - ik) 



(3.51) 



(3.52) 



Finally, inserting (3.51) into (3.50) and taking the limit n — > — oo we derive the 
evolution of spectral data. Indeed, evaluating ( [3.50 ) on the real axis and taking into 
account ( f3.14 ) and (2.16)— ( 2.19| ) we have 

n+l 



p ( 1 — ik 
2 [l + ik 



Therefore 



b+ ( k ) i — L n b+ ( k ) + C 1 + lk ) a t( k ) 



{l+ik) 



l + ik 



b+(k) 



(l-ik) n +\ +n ^ 

q-T, 777 70 k. 

(1 +ik) n - 1 



a+(t;k) = a + {0;k) 
b+(t;k) = 6+(0;fc)exp 



ik P-M\+k*) t 



Analogously one gets for the evolution of the normalization coefficients C + 



. f+ p-2?(l 
ikj 



k ? 



Cf(t) = C+(0)exp 
If p = 2q the evolution for b + (t; k) and Cj~(t) reduce to 
b + (t;k) = 6+(0;fc)exp 



k 3 



C+(i) = C/(0)exp 



'.+3 



-2ig 



1 



kf 



(3.53) 
(3.54) 

(3.55) 

(3.56) 
(3.57) 



which is consistent with the expected time evolution of KdV in the continuous limit. 
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4. Difference-difference equation 



We discretize now also the time and we consider the following Lax pair 

■020 = A?/>oo + 500^10 (4.1) 
ipu = (A 00 + XB 0Q ) ipoo + (Cio + AL>io) Vio- (4.2) 

If we denote by E and F the shift operators Efy = fi+\j and Ffy — /y+i and 
impose to (4.1) and ( [1.2] ) the compatibility condition Eipn = -Fi/^o we get 

{A w + XB 10 ) ip w + (C 20 + A£>2o) [AV'oo + ffooV'io] 

= (-4-io + A-B_ 10 ) [-010 - ff-ioV'oo] + A (C 00 + A-D 00 ) V'oo 

+ 5oi [(-^oo + AS 00 ) ipoo + (Cio + ADio) V'io] 

and, then, equating the coefficients of ipoo, V'io and of the different powers of A we get 

D 20 = Aio (4.3) 

G20 — Coo = goiBoo — g-io-B-io (4-4) 

3oiA)o = 9-wA- W (4.5) 

-Bio - £_io = ffoi-Dio ~ .gooAm (4.6) 

-A w + A_io = 300C20 - ffoiCio- (4.7) 



We choose 



Du = 



(4.8) 

Bij = b (4.9) 

and, then, this system of nonlinear recursion relations can be explicitly solved. It is 
convenient to introduce the following notations 
#01 5-11 



Too 



- 1 



500 5-10 
Goo = (501 — 5-io) 



Po 



1 



5- 205-10 



1 



5-205-n V 5ii5oi 
We get for the coefficients of the auxiliary spectral problem 

tooAio + ba Q0 = 

Too 7oo Coo + bao o- o = 6t oG o 



where 



too = (5oo) 2 /?2o7io - (5oi) 2 /3io720 
°oo = 500710G20 — 501720G10 

J 5oi_ 

5ii 



c*oo 



-(1 + 700) 



5-io5-20 

and the nonlinear double discrete equation for g 



(4.10) 
(4.11) 
(4-12) 

(4-13) 
(4-14) 

(4-15) 
(4.16) 

(4-17) 



5oo5oi 



(5-io) 2 /?io7oo - (5-n) 2 Ara7io 



G 



21! 



+ (5oi) 2 (5-11) A)o720 - (5-10/501) (500) /32o7oo 



5-105-11 



(5oo) 2 /32o7io ~ (5oi) 2 /3io720 



Go 



G10 
= 0. 



(4.18) 
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If we let 

g(n,m) = g n (mh) 
t = mh, 

taking into account that in the limit h 

7(n,m) {g n g n -x)t 



h 

(3(n, m) - 
G(n, m) 



g n g n 
1 



g n -ig n -2 



(4.19) 
(4.20) 

(4.21) 



(4.22) 
(4.23) 



9n+ig n 
9n - g n - 

it is easy to check that this equation in the limit of continuous time is satisfied by the 
differential-difference equation (1.2) which contains two arbitrary constants p and q. 
We expect that also in the general case the equation (4.18) could be integrated to a 
lower order equation containing two arbitrary constants. 
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